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Abstract 



Breathers are spatially localized and time periodic solutions of extended Hamiltonian dy- 
namical systems. In this paper we study excitation thresholds for (nonlinearly dynamically 
stable) ground state breather or standing wave solutions for networks of coupled nonlinear 
oscillators and wave equations of nonlinear Schrodinger (NLS) type. Excitation thresholds 
are rigorously characterized by variational methods. The excitation threshold is related to 



Hamiltonian energy. We establish a precise connection among d, the dimensionality of the 
■ lattice, 2a + 1, the degree of the nonlinearity and the existence of an excitation threshold 

Q\ ! for discrete nonlinear Schrodinger systems (DNLS). We prove that if a > |, then ground 

state standing waves exist if and only if the total power is larger than some strictly positive 
§ ■ threshold, v t hresh{o~,d). This proves a conjecture of Flach, Kaldko & MacKay |l3j in the 

context of DNLS. We also discuss upper and lower bounds for excitation thresholds for 
ground states of coupled systems of NLS equations, which arise in the modeling of pulse 
propagation in coupled arrays of optical fibers. 

> : 

1. Introduction 

f?, 

This article concerns threshold behavior of certain time- reversible, energy preserving nonlinear 

dynamical systems. Consider an infinite dimensional Hamiltonian system (wave equation or 

network of discrete oscillators) denned on an infinite spatial domain. If the system is translation 

invariant ( e.g., not having any localized potential wells), one expects that "small amplitude" 

or "low energy" solutions will disperse to zero; see, for example, ||24j| . If the system is nonlinear 

and having an attractive nonlinear potential, one can expect that sufficiently large "amplitude" 

initial data will lead to an evolution consisting of a non-decaying "bound state" component 
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plus a dispersive component (radiation), which tends (weakly) to zero with increasing time. 
In this latter scenario, we think of permanent non-decaying structures as having been excited 
by the initial condition; a deep enough self-consistent potential well has been initialized in 
which one can sustain a permanent structure. Since the systems we are discussing are infinite 
dimensional, the sense in which one measures amplitude is crucial. In systems of physical 
interest, there is often a natural measure of a solution's size. Roughly speaking, if there is a 
critical size, Vthresh > 0, such that there are permanent (non-decaying in time) states of size v 
if and only if v > v thresh , then we refer to v thresh as an excitation threshold. In this paper, we 
investigate the existence and nonexistence of excitation thresholds for a class of time-periodic 
and spatially localized standing wave states for two classes of dynamical systems. In certain 
models, these states have been called "breathers". See section 3 for a precise definition of and 
discussion concerning excitation thresholds. The dynamical systems we consider are: (1) the 
discrete nonlinear Schrodinger equation ( |1 . 1|) and (2) a system of coupled nonlinear Schrodinger 
equations (|1.7|) ; see also ( |6.1|) . 

Mathematical models which support discrete breathers are of interest in the study of vibra- 
tions in, for example, localized crystals and biological molecules 0, [14]. Recently, experimental 
observations of such discrete nonlinear localized modes have been made in coupled systems of 



optical waveguides [1C[. With a view toward study of such structures in experiment it is of 
interest to understand under what circumstances a discrete breather is excited. 



In [23] a formal variational argument is given suggesting the existence of such energy thresh- 
olds for the one-dimensional discrete nonlinear Schrodinger (DNLS) equation (also known as the 
discrete self-trapping equation M). For the related system of nearest neighbor coupled nonlinear 
Schrodinger equations, (|1.7|), such thresholds were rigorously demonstrated to exist [P^[ . 

In the recent paper of Flach, Kladko & MacKay fl3T| , heuristic scaling arguments and numer- 
ical studies are presented which suggest that for a large class of Hamiltonian dynamical systems 
defined on one, two and three dimensional lattices, there is a lower bound on the energy of a 
breather if the lattice dimension is greater than or equal to a certain critical value. 

Theorem 3.1 resolves this conjecture for ground state breathers of the d- dimensional discrete 
nonlinear Schrodinger equation (DNLS): 

idtipi = -k ( 5 2 if ) i - |^| 2<J fr, k > (1.1) 

Here, ip = {tpi(t)}, I G 7L d , t G IR, 5 2 denotes the rf-dimensional discrete Laplacian on 7L d 
given by: 

(<sV), = Ei - 2d ^ (i-2) 
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where iVj denotes the set of 2d nearest neighbors of the point in 7L d with label I. The parameter 
k can be interpreted as a discretization parameter, k ~ /?T 2 , where h is the lattice spacing and 
ipi = ip{hl), I G 7L d . The parameter a > is a measure of the degree of nonlinearity. 

Theorem 3.1 states that there is a ground state I 2 - excitation threshold if and only if a > |. 
For a < | breathers of arbitrarily small I 2 norm exist. See [|TTJ, \T2\, [I3~| for a study of the 



bifurcation of small amplitude states from the edge of the plane wave spectrum. In contrast, 
the continuum limit nonlinear Schrodinger equation, ( p.lp , has an L 2 threshold only in the case 
of critical nonlinearity, a — |. This is a manifestation of the role of discreteness, which breaks 
the dilation invariance of the continuum case; see the discussion and analysis of sections 3 and 
4. Theorem 2.2 states that ground states are nonlinearly dynamically stable in an orbital sense; 
see also |0J. 

In section 5 we consider the limiting behavior of ground states of total power. \\ip\\ 2 2 = v, 
as v tends to infinity. Such ground states are found to have large amplitude. As v is increased 
they are increasingly concentrated about one lattice site. A phenomenon of this type has been 
observed for the systems ( |1 . 7|) , (|6.1|) , and analytically studied in [27], The relation of this 



result to the numerical work of 0] and to the work on the anti-integrable limit [21 1, || is also 
discussed. 

Studies of discrete breathers originated in the context of classical nonlinear wave equations. 
An example is the one-dimensional Klein-Gordon equation: 

d 2 u n = D ( u n+ i - 2 u n + w n -i ) -Oj m„ + u 3 n . (1.3) 

The techniques of this paper do not directly apply to give rigorous thresholds for discrete 
nonlinear Klein Gordon equation localized states. However, our results concerning DNLS are 
related, through a multiple scale approximation, appropriate to the limit of large lattice spacing, 
h. Specifically, let h = e^ 1 , and therefore D = e 2 k. Then, seeking a solution of the 
form: 

u n = e ^> n + e 2 $ n + e 3 X n + ••• (1-4) 
we find an approximate solution which is valid for times, t, of order e~ 2 with 

# n = tf„(*,T) = e-^i; n (T) +e in ° t MT), T = h 2 t (1.5) 

where i/j n (T) satisfies the discrete nonlinear Schrodinger equation ( |1 . 1|) . In particular, this yields 
using the results of this paper on DNLS approximate solutions of the form: 

u £ n (t;u) = 2e cos ( [O + £ 2 u]t + 7 ) g n + 0{e 2 ), (1.6) 
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where u < and g = g u = { g n } neZ G 1 2 {7L). 

Finally, in section 6 we discuss and extend results on excitation thresholds for ground states 
of a class of coupled system of nonlinear Schrodinger equations (CNLS), which arises in the 
modeling of pulse propagation through a coupled network of optical fibers @, PD| , |p7f|, 
BOB: 



id t ^ + %<l>i + « ( S 2 if ) i + 2 |^| 2 ^ = 0, (1.7) 
$ = {ipi(t,x)} le%d , (t,x) G IR 2 

The cases of physical interest are d — 1,2. Here, ipi denotes the slowly varying envelope of 
the highly oscillatory electric field in the fiber with position I in the lattice. We consider the case 
where I varies over 7L d , with J2i H^IIz^r) < 00 • For the case & = d = 1, we obtain numerical 
upper and lower bounds (|6.23|) for the excitation thresholds z^ c .[] Other boundary conditions 



are discussed in ||27|| , |30| . In particular, a result of the analysis is that there are no excitation 
thresholds in the case when the system is periodic in the discrete variable, l; ground states of 
arbitrary positive total power v = Yli II ^11 1 exist . 

In this paper we use observations about the scaling structure of variational problems together 
with compactness methods in the calculus of variations; see e.g. f22|, p. Thresholds for the 



excitation of breathers or nonlinear bound states are characterized in terms of the optimal (best) 
constant of discrete interpolation inequalities for elements of l 2 {7L d ) in the case of (|1.1|) and for 
elements of I 2 (2L d ; H 1 (JR)) in the case of ( |1.7| ). This is related to the approach taken in [p8[j . 
29] on a sharp criterion on initial conditions for global existence (no finite time blow-up) of 



solutions to the continuum nonlinear Schrodinger equation on IR , ( |3.1| ), with critical power 
nonlinearity. Results on excitation thresholds, stability and other issues for the semi-discrete 



class of nonlinear Schrodinger equations were obtained by B. Yeary and the author ||27|| , [PCj] . 
This article is a detailed account with extensions of the work on excitation thresholds. 

Acknowledgements: The author wishes to thank Peter D. Miller for very helpful com- 
ments. This work was supported in part by a grant from the U.S. National Science Foundation. 

2. DNLS and a variational characterization of its ground state 

By standard methods, one can check that for any ip(t = 0) G l 2 {7L d ) ) there is a unique global 
solution if G C 1 (JR; I 2 CZ d )) of DNLS, (|1.1|) , and for which the following two quantities are 
^^An error in these bounds due to faulty algebra appeared in |2^] and is corrected here. 
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independent of time: 



TCd 



1> 



- K (s%^) — — Yl M 
v ' <j + 1 . 

E m 2 - 



2a+2 



(2.1) 
(2.2) 



2GZ d 



The subscript, " D" , is used to indicate a quantity associated with the discrete equation (1.1) 



TLd is a Hamiltonian for (|1.1|) , which can be written as: 

6H D 



id t ij) 



5i/j* 



(2.3) 



In various applications the invariant TV has the interpretation of total power or of particle number. 
The term — may be written out explicitly as: 



= E E hfo-ift* 

r=l iGZ d 



(2.4) 



where r r denotes translation by one lattice unit in the r" 1 coordinate direction. 

Of particular interest are spatially localized and time-periodic solutions. We seek them in 
the form: 



(2.5) 



where is real. A solution of this type is frequently called a nonlinear bound state, standing 
wave or stationary state. The term discrete breather is also used but is sometimes reserved for 
a localized state whose modulus oscillates. 

Substitution of (|2.5|) into Ql.lD yields the system of algebraic equations plus "boundary 
condition at infinity": 



w 9i 
9 



(2.6) 
(2.7) 



We construct a ground state by variational methods. To motivate our approach, we consider 
the quantum mechanical problem: 

H * = E ip, (2.8) 
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where H = — A+V(x) for a bound state, \1> G L 2 with H^l^ = 1- We assume V(x) is a sufficiently 
smooth and rapidly decaying "potential well". Consider the constrained minimization problem: 

J = inf{ (HfJ) : ||/|| 2 = 1 }. (2.9) 

If X < 0, then E g = 1 is the ground state (lowest) eigenvalue and there exists a ground state 
eigenstate ^f g (x) such that 

Ht> g = E g V g , ||ttj 2 = l. (2.10) 

The time-periodic breather or standing wave, \& g (x)e~ lEat , is a dynamically stable solution of 
the time-dependent Schrodinger equation 

id t V = H (2.11) 

We shall characterize the ground state of ( |1 . 1|) using a nonlinear analogue of ( |2.9|) . 

Definition: Let 

J, = inf{H D [ /] : AT D [ f] = u}. (2.12) 
A minimizer of the variational problem Q2.12 ) is called a ground state. 



Clearly, X v is bounded below: For 

I 



n D [f] >-— rEi^r +2 >--— rll/llL CT H/ll2>--— r^ +1 - (2-13) 

(7+1, ct + 1 cr + 1 



Theorem 2.1. (a) If — oo < Z„ < 0, then the minimum in ( |2.12| ) is attained. 



(h) Every minimizing sequence associated with the variational problem ( |2. 12\ ) is precompact 
modulo phase translations, i.e. for any minimizing sequence { (f® }, there is a subsequence 
{ cf n ^ } and a sequence { 7 nfc }, and translations, r(l k ) (where r(l k )g^ = { g^i} j e z d )> sucn 
that r(/fc)^" fe ) e lJn k converges in l 2 {TL d ) to a minimizer. 

(c) If g = {g/} igZ d is a minimizer for the variational problem ( 2.1% ), then there exists u> = 
oj{v) < such that the Euler-Lagrange equation: 

u(u) 9l = -k (S 2 g) i - \gi\ 2a g h le7L d (2.14) 

holds, together with the L 2 constraint: 

Md[9] = El^l 2 = (2-15) 
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This theorem can be proved by a standard application of concentration compactness ideas 
in the discrete context [^2|; see |30| . An outline of the proof is presented in appendix A. 



Dynamical stability 

Before stating a precise result, we first introduce some terminology and notation. 
Definitions: 

(1) Let Q v denote the set of all solutions of the minimization from ( 2.12p , i.e. the set of ground 
states with J\f — v. 

(2) Given a particular ground state g, we define its orbit to be the set: 

0{g) = {^g : 7 e [0,2tt)} (2.16) 

(3) The distance p ( ip , Q v ) from tp e I 2 to the set of ground states, Q v is given by: 

p($,Gu) = jnf || f - <? |U Z d) (2.17) 

Remark: We conjecture that the ground state with M = v is essentially unique, i.e. if g is any 
ground state with M\ g] — v, then Q v = O(g). 

A consequence of part (b) of Theorem 2.1 is the following J7|: 



Theorem 2.2. Ground states of ( |1 . 1\ ) are orbitally Lyapunov stable in the sense that: given 



any e > 0, there is a 5 > such that if the initial data ip(t = 0) — ipo satisfies 

p($o,G„) < 5, (2.18) 

then for all t ^ 

p(m , Q v ) <e. (2.19) 
3. Excitation thresholds for DNLS 



For a fixed lattice dimension, d, we consider the family of equations (|1 . 1|) parametrized by a. 
Theorem 2.1 gives a criterion for the existence of a ground state. 



Definition: If for any v > the variational problem ( |2.12| ) has a strictly negative infimum 



X v < then, by Theorem 2.1, a ground state exists for any v > 0. In this case we say that there 
is no excitation threshold. However, if there is a strictly positive constant y^ res h (which may 
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depend on d and a) such that X v < if and only if v > ^ res h, then we call vf hresh an excitation 
threshold or L 2 excitation threshold for a ground state. 

The main result concerning DNLS is the following: 

Theorem 3.1. (1) Let < a < |. Then, I u < for all v > 0. Therefore, the variational 
problem ( 2.12 ) has a solution for all v > and there is no excitation threshold. 



(2) Let a > |. Then, there exists a ground state excitation threshold, v^ Tf , sh > 0. 

Remark on DNLS vs. NLS: Here we contrast the discrete equation, DNLS, and its continuum 
limit. In particular, we comment on some consequences of the breaking of various symmetries 
in passing from NLS to DNLS. 

(1) The continuum limit of ( |1.1|) (k = h~ 2 , h = lattice spacing, and h — > 0) is the d— dimensional 
nonlinear Schrodinger equation: 

id t <j) = -A (p - |0| 2a 0, (3.1) 

For initial data <f>(t = 0, x) G i^ 1 (IR d ), it has been shown that there exists a local solution which 
is continuous in time with values in H x (IR d ) and which satisfies the analogous conservation laws 



| 16|j , ||18|| . If a < 2 solutions are always global in time, while for a > |, finite energy initial data 



may give rise to a solution which leaves the space H\]R d ) after a finite time 0, @- In 
contrast, the evolution for (|1 . 1|) is globally defined in time. 



(2) Solitary standing waves can be found by methods analogous to those used in section 2. An 
excitation threshold for standing waves, in terms of the natural L 2 invariant: 

= I J \m\ 2 dx (3.2) 



exists only in the case a — |. This follows because under the scaling: 

(j)(x,t) ^4> p (x,t) = p-(/){px,p 2 t), (3.3) 

we find 

M[<f> p } = pi-'Micj)}. (3.4) 

Thus, given that a single standing wave exists, if a ^ |, scaling can be used to find one of 
arbitrarily small total power, Mnls- In contrast, the dilation symmetry is broken in the discrete 
case. 
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(3) Let a = |, and let R denote the ground state standing wave. That is, R is an H 1 solution 
of A R — R + Rd +l = of minimal power Nnls = Nthresh- In 0], ft was proved that if 

Af[0 o ] < Mthresh (3-5) 

then the solution exists for all time and disperses to zero in the sense that ||0(t)||i> — > 0, as 
|t| — > oo, for p > 2. 

Conjecture: 7/A/"[ V'o ] < ^thresh' then the solution of DNLS disperses to zero in the sense 
that for any p G (2, oo] : 

||^)L(z d ) 0, as |t| -> oo. (3.6) 



Proofs of the assertions in (3) are given in p9f,|28[ and rely on the pseudo-conformal symmetry 



of the continuum limit NLS, a symmetry which is absent in DNLS. 

Theorem 3.1 is a consequence of Propositions 4.1 and 4.2 of section 4. We begin by investi- 
gating the conditions on a, d, and v under which X v < 0. 

Proposition 3.1. I u > if and only if is is such that the following inequality holds for all 
uel 2 {7L d ): 



E N 2CT+2 < + MT H 2 (-6 2 u,u). (3.7) 

To prove Proposition 3.1, we observe that X v > if and only if for all u G l 2 (7L d ), with \\u\\ 2 2 = v 

(a+1)- 1 E \ Ul \ 2 °+ 2 < n(-5 2 u,u). (3.8) 

IC7J' 

Let 7^ v G I 2 be arbitrary. Then, if u defined by: 

u = y/u \\v Wp 1 V (3.9) 



satisfies the inequality (|3.8| ), which after some algebra yields (|3.7|) . Finally, if > we 
have that X v = 0. This is seen by simply taking a sequence whose N th element is a constant 
(depending on v) on the set of sites satisfying |Z| < N and zero otherwise. Along such a sequence 
we have M — v and TC tending to zero. Therefore, I u = 0. 



9 



Strategy of the proof of Theorem 3.1: Clearly, if the inequality ( |3.7|) holds for some v\ 
then it holds for all v < V\. We shall prove in Proposition 4. 2d that a ground state does not 
exist for any < v < v\. We are interested in characterizing fthresh defined by: 

^thresh = sup{ v : inequality (3.7) holds}. (3.10) 

In the following section we relate this threshold value to the optimal (best) constant in an 
interpolation estimate related to the Hamiltonian energy. H. If a finite positive vf hresh exists, 
then for any v > v^ hresh and element of l 2 {7L d ) ) w*, can be found which violates the inequality 
( |3.7|) . This choice of u* shows that X v < 0, and by Theorem 2.1 there is a ground state. If, 
however, for any choice of v > one can construct an element of I 2 for which the inequality (|3.7| ) 
is violated, Theorems 3.3 and 2.1 imply that a ground state exists for any v > 0, i.e. there is no 
excitation threshold. The strategy used to prove Theorem 3.1 is to show that if < a < |, 
then there is no value of v for which the inequality ( |3.7p holds for arbitrary u e I 2 . However, 
if a > | we show it holds if and only if v < v® hresh , for some v® hresh > 0. 



4. Best constants and excitation thresholds for DNLS 

In this section we relate the problem of characterizing excitation thresholds to the problem of 
finding the optimal or best constant in discrete interpolation inequalities of Sobolev-Nirenberg- 
Gagliardo type. 

The discussion concluding section 3 motivates the following question, answered in Theorem 
4.1 below: 

When does there exist a constant C > such that for all u = {u{\ E l 2 {7L d ): 

E KI 2CT+2 < C ( £ \uA (-5 2 u,u)? (4.1) 
h}7J l \ h:7J l J 

If ( [4.1|) holds for some C > and C* is the infimum over all such constants, then v® hresh 
defined by 

(£7+1) « (ug resh )~ a EE C. (4.2) 

is a ground state excitation threshold. Therefore, we seek to characterize the optimal constant, 
C*. If there is a strictly positive and finite C, then 

i = j<t,d s M (£ ^; |M ' |2 >"(;f""""> (4.3) 
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and we have: 



^thresh 



( (a + 1) k J a " 



(4.4) 



Remark: If J a,d > 0, then by Proposition 4.2 below, there exists a strictly positive lower bound 
on the energy, Af, of a ground state. 

Note that ( [4.4Q is consistent with the simple observation that for the case of uncoupled lattice 
sites, k = 0, there is no excitation threshold. For example, in this case the solution 



Mt) = 0/^0 (4.5) 
is a l 2 {7L d ) solution of (|1 . 1| ) with M = v. This limit is also called the anti-integrable limit |^l] 



||]. In section 5 we shall relate the anti-integrable limit to the large amplitude limit of our 
variationally constructed ground states. 

Proposition 4.1. If a < §, then J a,d = 0. Therefore, for a < §, and there is no ground state 
excitation threshold (v 'thresh = ®)- ^ D other words, ground states of arbitrary energy, M , exist. 

proof of Proposition 4-1- Consider the one parameter family of trial functions, u(a) defined by: 

u t (a) = e~ a I'l, (4.6) 

where / = (h, ...,ld) £ ^> d , \l\ = |^i| + ••• + \ld\ an d a > 0. Evaluation of the terms of the 
quotient in ( |4.3|) yields, for a | 0: 



2 -d 

~ a 



a 



2-d 



7 K 



2(7+2 „ -d 

~ a 



(4.7) 



Therefore, the quotient in 
This proves the Proposition 4.1. 

Proposition 4.2. Let o > ~. Then, 

(a; > o. 

(b) If || ^ ||f a = i/ ; then 



is of order a 2 dcr , which tends to zero as a tends to zero if a < |. 



v 



thresh , 



(4i 



where v^resh > JS given by (fOp. For a > |, ^threshi '^ d) JS azi excitation threshold, i.e. 

(c) if is > vf hresh {p, d) then X u < and a ground state exists, and 

(d) if v < v^hreshk ^ c 0; *-^ eii = and there is no ground state minimizer of ( |2.12| ). 
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proof of Proposition 4-2 : To prove part (a) suffices to show that the inequality ( |4.1| ) holds for 
some positive constant, C. Then part (c) follows from the discussion at the end of section 3. 
We proceed as follows. For functions / G i^ x (IR n ), one has the Sobolev-Nirenberg-Gagliardo 
inequality |0J: 

mini < c iiv/ni/ii2 +<j(2 " n) , (4.9) 

where a is restricted to satisfy: 



< a < oo, n = 1, 2 

< a < 2(n~2y\ n > 3. 



(4.10) 



The proof of ( f4.9| ) can be followed closely to yield, under the same restrictions on a, the following 
estimate in the discrete case for u G l 2 {TL d ): 

l+f (d-2) 



E M 2ff+2 < c £ \ Ul \ 



( —5 2 u, u ) 



ad 
2 



(4.11) 



To give the idea, we present the proof of (|4.11|) in the case d = 2. We write u\ = u ab , (a, b) G 71? . 
Without loss of generality we can take u ab > 0. Note that 



u 



.<T+1 

ah 



E 



u 



.(7+1 

ab 



II 



cr+1 
a-l.b 



(4.12) 



By the fundamental theorem of calculus 

-i d 



u 



(7+1 

ab 
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(T+l 
a-l.b 



/ T~ [ SU <*b + (1 - S)U, 

Jo ds 



a— 1,6 



icr+l 



ds 



o- 



1) 



SU ab + (1 - S)u a _ hb Y ds ( U ab - lta-1,6 ) 



Therefore, (using the convention that sums without specified upper and lower limits are under- 
stood to be taken over all 7L d } 



.(T+l 

ab 



U 



.cr+1 
a-l,b 



< \(T+1\ E max(|w a6 r > l^a-l^D I «a6 - Ma-1,6 I 

It follows by summing over a and applying the Cauchy-Schwarz inequality, that 



E 



u ab u a-l,b 



<2^|a + l| Kb\ 



2a 



(4.13) 



(4.14) 



The analogous computation can be performed by summing on the second index, to get: 



E US 1 ^ <J-i < 2*k + i| E E K P - 



u a,P-l\ 



(4.15) 
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The product of the last two estimates yields: 

< 2|a + i| 2 rr m 2ct 



\Uab 



2a+2 



E \ u ab — Ua-\,b\ 



X 



E K, 



|2(T 



E \u a p- Ua,{3-1 



Summing on a and applying the Cauchy-Schwarz inequality gives: 



E k&r +2 < 21^ + 1| 2 hr k & M e i«o6 



E K 



■01 



|2(T 



E K/3 ~ Wo,/3-l| 



Finally, summing this result on b and applying the Cauchy-Schwarz inequality gives ( 4.11[ ) for 
the case d = 2 and arbitrary a > 0. 

To complete the proof of Proposition 4.2, we write estimate fl4.11|) as: 



E M 2<T+2 < C E (S 2 n,n 



6 u, u 



ad -i 

2 — 1 



IC7J' 



,J2lcZ d \ u l\ 



(4.16) 



The last factor in ( |4.16| ) is bounded by a constant for a > |; the discrete Laplacian is a 
bounded operator. Therefore, if in addition to ( 4.10|) , we have a > |, then the estimate ( |4.1| ) 
holds. 

Finally, we want to show that for d > 3 we can relax the constraint < 2 — a(d — 2). 

Suppose d > 3 and take a in the range for which we know the estimate Q4.1| ) to hold. This 
estimate is equivalent to: 

E N 2ff+2 < C (S 2 u,u) (4.17) 



subject where u satisfies the constraint 



E M 



1. 



(4.18) 



The constraint ( 4.18Q implies that for all I G 7L d , \ui\ < 1 and therefore if a±, is any number 
satisfying o\ > a > §, then the estimate (^4.17|) holds with o replaced by o\. This implies the 
following result which completes the proof of part (a) Proposition 4.2: 



Theorem 4.1. For a > |, the interpolation inequality holds. 
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Remark: Note that there is no upper restriction for n > 3 on a as in the continuum case 
Through ( |4.16| ), the boundedess of the discrete Laplacian, — 5 2 , on l 2 {7L d ) plays a key role. 

Part (b) of Proposition 4.2 follows from the definition of Tip and the inequality (|3.7Q with 
optimal choice v = ug resh given by 

Finally, we prove part (d) of Proposition 4.2. Suppose v < ^ resh - Then, by part (b), X v > 0. 
On the other hand, as at the end of the proof of Proposition 3.1, we have that X v < 0. It follows 
that I u = for any v < v® hTesh . If the minimum is attained at a state tp, then 



a + 1 l 

E m 2 = " 



Since cr > |, v® hresh defined by ( |4.4[ ) is strictly positive and by (p.7[), with the optimal choice 

v = "thresh, w e have 

k ( -5V , ^ ) < « f"^— 1 ( -<^V > ^ ) < « ( (4-19) 

V U thresh ) 

a contradiction. 

Theorem 3.1 now follows from Propositions 4.1 and 4.2. 

5. Large amplitude and the anti-integrable limit 

Large v limit of ground states 
As discussed in [[21]], breather solutions of DNLS can also be constructed perturbatively in 



the limit of zero coupling, k = 0, also called the anti-integrable limit; see also ||. In [f2~l"l 



as an explanation for the numerical studies in |4[], it is conjectured that the large amplitude 
anti-integrable limit breathers play an important role in the dynamics of DNLS. We now give 
evidence of this, by showing the connection between the nonlinearly stable ground state breathers 
constructed by variational methods and the large amplitude anti-integrable breathers. We also 
prove that as v increases, ground state breathers of "total power" v grow in amplitude and 
become increasingly concentrated about one lattice site. This property of ground states and 
their nonlinear stability (Theorem 2.2) elucidate the numerical simulations in 

We begin by considering a scaled version of the variational problem ( |2.12| ), for the DNLS 
ground state: 

^ = inf{ -K(5 2 f, f) - -L- £ \fr +2 ■ E \fi\ 2 = " }• (5-1) 
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In anticipation of our taking v f oo, we set 

/, = V \ F h £ W = !• (5-2) 
i 

and introduce the parameter 

= -, 5.3 



which tends to zero as v tends to infinity. The variational problem, (|5.1|) , is then equivalent to: 

K{a) = inf{ a (-5 2 F, P) ■ -J- E |F Z | 2 ^ 2 : E |F,| 2 = 1 }. (5.4) 
V J a+1 i i 

By Theorem 3.1, if v > vf hreah > 0, then there is a ground state breather: 

G = G(a) = {G,(a)} l6r i, (5.5) 
satisfying the Euler-Lagrange equation: 

-a(5 2 G) i -\G l \^G l = \ a G h (5.6) 



where X a is a Lagrange multiplier. By ( |5.2| ), this gives rise to a ground state family of solutions 
of O): 

f g {t) = i/s G(a) e - 4A ^ CT * e* 7 , 7 G [0,2tt). (5.7) 

By Theorem 2.2, the ground state family is nonlinearly orbitally stable, and is therefore expected 
to participate in the dynamics. 

What is the structure of ground states for large v? We next show that as v — > oo ; ground 
states become concentrated on the lattice about a single site. 

To see this, we first observe that by the methods of appendix A (see also |§), 
(a) As a tends to zero {y \ oo) through a sequence, {G(a)}, is a minimizing sequence for the 
limit variational problem: 

1 



inf{ -— E \Gi\ 2a+2 ■ E l^l 2 = 1 }• (5i 



(b) A subsequence can be extracted, which (modulo phase adjustments) converges to a mini- 
mizer, G°°. G°° satisfies the Euler-Lagrange equation associated with the limit problem ( |5.8|) : 

GOO I 2(7 /OOO \ /~ioc 

, I \ u i — A , 

E \GT? = I- (5-9) 
I 
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7,1 r ' x - " ll • ' " • - - r " -> 1 Since || G°° || F(z<i) = 1, G 1 * 

can be nonzero only at a finite number of sites, N > 1. Therefore, 



Thus, for each / G G z °° G {0} U {(-A)*? e* 7 : 7 G [0, 2n) }. Since || G°° ||, a(B<l) = 1, G^ 



|| G°° || ;2(zd) = 1 implies - A = A^ ff 

1 |£-foo|2cr+2 iv 

(7 + 1 1 (7+1' 

The minimum is therefore attained for N — 1 and we have: Gf° = ±5// for some Iq G and 
A = -l. 

Therefore, as v — > 00, a subsequence of ground states converges, to a limiting state: 

G,(a) - G?° = 5 hl0 , in / 2 (S d ) (5.10) 

for some Iq G Therefore the large v [a small) limit of ground states behaves as a large 
amplitude one-site breather. 

frit) ~ ±1/* 6 l>l0 e~ iv °\ for !/ large. (5.11) 

Connection with the anti-integrable limit 
For small a equation (|5.6j ) is the anti-integrable limit studied in |^TJ. The approach taken in 



21] , |3j is to first observe that for a = each lattice site evolves independently and that ( |5.6| ) 
has solutions ^i(t), where for each I G 7L d , satisfies the equation: 

idtV(t) = -|^(t)| 2CT *(*). (5.12) 

The solutions of (|5.12|) are: 



= u e* M 1 e'\ (5.13) 

with w,7 G E. Fix a solution which is supported at lattice sites q G / C 7L d , where I is finite 
or infinite and such that at each site the evolution is an oscillation of the form 

tf ff (t) = u q e iM2at qel (5.14) 

and such that the frequencies uj q are all commensurate. The implicit function theorem implies 
that these solutions have a continuation for a sufficiently small in the space of time-periodic 
solutions. These range in spatial complexity from those that are small perturbations of the 
simplest a = breather, consisting of a solution of the form (|4.5| ), to those which are small 
perturbations of an a = "spatially chaotic" configuration of oscillators. 



16 



Consider the continuation from the anti-integrable limit of one-site breathers. These are 
solutions of the form: 

%(t,a,fi) = Ma^)e- ip \ I E TL d (5.15) 

where 

A l (a = 0,/j) = (-n)&8 ltlo , for some l E 7L d . (5.16) 

We wish to relate the two-parameter family A{a, /i) to the family of scaled ground states G(a), 
for small a. Note that A{0,—1) is such that ||A(0, — l)||ja(z d ) = 1- It is easy to check, by the 
implicit function theorem that a locally unique solution (a,/i(a)) defined in a neighborhood of 
a = exists such that 

\\A<*> K<x))\\p(x*) = h 

MO) = -i. 

Therefore, by our variational arguments and local uniqueness: 

G(a) = A(a,fi(a)). (5.17) 

6. Thresholds for coupled systems of nonlinear Schrodinger equations 

In this section we discuss results for systems of coupled on nonlinear Schrodinger equations 
(CNLS) (0): 

i) = {^/(t,x)} ieZd , d=l,2, (t,x) G IR 2 

CNLS has been introduced as a model governing the propagation of light pulses in a coupled 
d = 1 or d = 2 dimensional array of optical fibers. We consider the case where the discrete 
variable varies over 2L d , and such that ipi(t, x) decays as I and x tend to infinity. Other boundary 
conditions {e.g. periodic are considered in [ppf] , §J. We follow a similar outline for the CNLS 
as that followed in our analysis of DNLS. Certain details are omitted and for them we refer to 

Given initial data ipo{ x ) f° r CNLS satisfying 

Y, < oo, (6-2) 

lrZ d 
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there is a unique solution t i— > ip(t,x) which is continuous in t with values in l 2 {7L d ) x i^ x (IR). 
The following two functionals, evaluated on solutions, are independent in time: 



H[$] = I (-5 2 j(x)J(x)) + J2 I \dM*)?fc- I \ip(x)\ 2a+2 dx 
= E/ \in?dx. 

TC is a Hamiltonian energy of the CNLS in the sense that CNLS can be expressed as: 

sn 



id t tp 



(6.3) 



The functional M corresponds to the total input power in the system. 

Of interest are nonlinear bound states of CNLS. These are solutions of the form: 



if = e iA * g(x; A) 



(6.4) 



for which the invariants Ti and M are finite. The components of ip satisfy the coupled system 
of equations: 



-A 2 ^ + 8% + k (5 2 $) l + (a + l)|^| 2<T ^ = 0/ 



G 7L° 



(6.5) 



In analogy with the discrete case, we seek to characterize the ground state of the system by 
variational methods. 
Definition: Let 



J v = inf{H f 



n 



f 



(6.6) 



Because of the similarity of the arguments to those in the previous sections and the more detailed 
treatment in p7|,[p0l we provide a summary. 



(1) If < a < 2, then J v > -oo for any v > 0. 

(2) In analogy with Theorem 2.1, we can show: 



Theorem 6.1. The infimum in 6 ) is attained if and only if J v < 0. Moreover, any minimizing 
sequence has a subsequence which converges strongly in I 2 {TL d \ H 1 (TR)^J modulo translations in 
space and phase. Furthermore, any minimizer satisfies the equation flO[ ) . 



In view of this result, we study the question: for which a, d and v do we have J v < 0? 



(3) 
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Proposition 6.1. J v = if and only if for any if) G I 2 \ 7L ; H l {Wi)j we have the estimate: 

II $ K2 < V ~°<£~ X II $ Iff <-*V . ^ ) 1_f II ^ II2, (6-7) 
where a a = ff ) 2 ~ CT — (f ) 



2 

2-ct 



1 



Here, ll/lg = ( E< II /< IIlp(r) j P - 

Proposition 6.1 is proved by a simple scaling argument. For any -0, such that || ^ III = v i 
we define the scaling ip r (x) = r^ip(rx), which preserves the L 2 norm. Evaluation of the 
Hamiltonian on ijf and minimization over r > gives: 

UW} >HW™\ = (-5 2 j, $) - a a || \\dJh^ (6.8) 

We can pass to an expression for arbitrary ijj by replacing ip by 1/2 1| ip W^ 1 ifj in (|6.8|) . This gives 
for any if G / 2 (V; F^K)): 

2 + CT 

W W ] = ( ^jj) <-*V, ^) - a, (^J II ^ ll^lMla^ (6-9) 



It follows that J v can be realized as the infinimum of the expression in ( |6.9| ) over all ip G 
I 2 (7L d ; H^JR)). Thus J v > if and only if (0) holds for all ip E I 2 (ZL d ; H\M)). As in the 
discrete case, it is simple to construct a sequence along which the L 2 constraint is satisfied and 
the Hamiltonian tends to zero. 

(4) Suppose an estimate of the type (|6.7| ) holds. In particular, we let denote the smallest 



constant for which this estimate holds. That is, 

II i> K2 < <?* || i> ||f (-5 2 if, f I) dj \\°, (6.10) 

where 

c , = ^ s . mf \\lR tfjjri 11 aJK (611) 

II Wj || 2o "+ 2 

II r II 2cr+2 

There are two possibilities. First, if K, (T ' d = C^ 1 = 0, then for any v > 0, there is a choice 

of ip which makes the Hamiltonian negative. In this case, by assertion (2), a ground state of 

any prescribed L 2 norm exists; there is no L 2 - excitation threshold. The second possibility is 

——l 

that < C7 1 = /C CT ' d < 00. In this case, we have that J v > if and only if C* < v~~ a aZ 
Therefore, we can define the threshold power, u c = u c (cr, d) by: 

u c = aV° C*° = aV° (/C CT ' d )". (6.12) 
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Use of the estimate ( |6.7| ) with the optimal choice v = v c , we obtain the sharp lower bound 
for the Hamiltonian, in analogy with the discrete case (compare (|4.8|) ): 



Kr\ > {-,Fr. (I (£) )• (0.13) 



for any tp with || tp \\ 



(5) The question of when an L 2 threshold exists is reduced to the determination of the range of 
values of a and d for which K. a,d > 0. Formula ( |6.12| ) then gives an expression for the threshold. 
To determine when JC cr,d is strictly positive amounts to determining when one can prove an 
inequality of type ( |6.10 ) for some (not necessarily optimal) choice of C*. This is addressed 



in 



[30f . Ranges of a, d for which this inequality fails to hold are determined by scaling 
arguments, while a proof of such inequalities for certain a, d can be obtained following the 
strategy used in the fully discrete case, where we mimic the proof of continuum interpolation 
estimates (e.g. see the proof of Proposition 4.2) or alternatively by applying the continuum 
interpolation estimates to functions of d + 1 variables and where the functions are taken to be 
piecewise linear in the variable corresponding to the d discrete variables; see f27|, [|30|. The 
results obtained are that K a ' d > for all a G [1,2) and for all er e (^ry, dZj)- I 11 summary we 
have: 

Theorem 6.2. Let a G [1,2) or a E (^tj, dZj)- Then, there exists an L 2 excitation threshold 



given by v c in ( \6.1%) . 



Estimates on v c for a = 1, and d = 1 
We now consider the case a = 1 and d — 1, an infinite one-dimensional array: 

id t ipn + d 2 x ip n + k( il) n -i - 2tfj n + tfj n+1 ) + 2 |^ n | 2 ijj n = 0, n e TL (6.14) 

We show how to get upper and lower estimates for the threshold power. A sketch was given 
in p7| , where an error appears in the displayed upper and lower bounds (due to an error in 
algebra) . 

By the above discussion, we know that there is an L 2 excitation threshold. That is, there is 
a constant v c > such that there are no ground states of with L 2 norm less than and there 
are ground states of L 2 norm z/2 for any v > v c . By ( |6.12|) (using that we must replace 5 2 by 
k,5 2 we have 

, e (l,l;«) = 2k\K* = 2/ci inf \\lM^LJ)llMh . (6 . 15 ) 

II ip lit 
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Upper estimate on p c (1,1;k): 

An upper estimate is obtained by evaluation of the functional in (|6.15|) on any ip ^ 0. In 
particular, if we use as a trial function the exact one-soliton supported on one site of the lattice, 
ipj(x) = sech(x)<5jo, we obtain the upper bound 

i/ e (l, 1; «) < 2v^6 ~ 4.89... (6.16) 



Lower estimate on vJl, 1; k): 



To obtain a lower bound we follow the strategy in |p7|| . First note that for arbitrary functions 



< Csng ( \\d x ip \\l + || d y ijj \\t ) || V HI- (6.17) 



By scaling in y, ip(x, y) i— > ^(x, ry), we have from (|6.17|) the estimate: 

II V lit < ZCsng II Ha II ^ lb II ^ 111, (6.18) 

In the best constant in ( J6.17|) is calculated and was found to be: 

C SNG = (tt x 1.86225.. (6.19) 

By ( |6.15| ), to obtain a lower bound for v c it sufficies to obtain an lower bound for JC 1 ' or 
equivalently an upper bound for C*. 

We next relate C* to Csng- This can be done by considering ( |6.17| ) for the restricted class 
of functions, tp{x,y), which are smooth in x and piecewise linear in y with jumps in d y i(j(x,y) 
at the integers. In particular, let 

ip(x,y) = (1 - 9) ip n (x) + 9ip n+1 (x), 
y = n + 0, < 9 < 1. 

Direct calculation gives: 

5 53 / l^n(^)| 4 ^ < y ^(x.y)^ dx dy 
Yl I \^n{x)\ 2 dx > / \if)(x,y)\ 2 dx dy 
/ \d x ip n {x)\ 2 dx > / \d x ifj{x,y)\ 2 dx dy 
Y. I I <Pn+i(x) - M*) ? ix = 

(6.20) 
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This, together with ( p.!8| ) yields: 

II $ lit < 5Csng \\f || 2 . (6.21) 

with a non-optimal constant, C = 5Csng which is an upper bound for C*. Thus, 

v c {1,1,;k) > 2 C; 1 > 2 {hCsNcY 1 

1 2 1 
= —Ti x 1.86225... > k* 2.3402... (6.22) 
5 



Combining (|6.16 ) and ( |6.22| ) we obtain: 



1 1 



2.34... < i/ c (l, 1, ; «) < «3 4.89... (6.23) 



A careful numerical simulation [W\ indicates k 2 v c (1,1;k) ~ 4.08. 



7. Appendix - Concentration Compactness Methods for DNLS 

Theorems 2.1 can be proved using the concentration compactness principle; see, for example, 



22]] . Since arguments follow quite closely those for the continuum case (see | 30fl , f22|j for a 
detailed implementation), we present here an outline of the ideas. 
Let = {ui}, denote a sequence in l 2 {TL d \ and such that 



EKT = v- (7-i) 
1 

Let B t (m) denote {/ G 7L d : \l — m\ < t}, and the norm \l — m\ = max!<j< d |/j — rrii\. 

Theorem 7.1. (Concentration Compactness Principle) 

There exists a subsequence satisfying one of the following three scenarios: 

(1) Compactness (the "mass" of the sequence concentrates): There exists G TL d such that 
for every e > 0, there exists a real positive number R e (independent of k), such that 

E K K) ! 2 > "-e (7-2) 

leB R£ (m k ) 

(2) Vanishing (the sequence spreads its mass over larger and larger sets and tends to zero): For 
all R < 00, 

lim sup J2 \ u t k) ? = (7.3) 
fc ^°° ™ez d ieB R (m) 
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(3) Dichotomy (the sequence concentrates its mass in at least two regions which become increas- 
ingly distant): There exists a G (0, v) such that , for all e > 0, there exist ko > 1 and disjointly 
supported sequences oS \ b (k) in l 2 {7L d ) satisfying for all k > k : 



\n^\\ 2 , 



U 2 (Z d ) 

— a 



2 

l 2 (Z d ) 



v - a 



< e 

< e 

< e 



distance (^supp(a ( - fc \ supp(6 ( - fc - ) 



OG 



as k 



oo. 



To prove this result, we introduce the sequence of concentration functions: 

Q^ k \t) = sup £ h (fc) | 2 , (7.4) 

meZ d l£B t (m) 



By following the arguments in P2| it can be shown that: 

(A) along a subsequence — > oo, Q nk (t) converges to a nondecreasing and nonnegative function, 
Q(t) with limit: 

lim Q(t) = a E (0,u). (7.5) 

t — >oo 

(2) the cases a = 0, a = v and < a < v correspond, respectively, to the above scenarios: 
vanishing, compactness and dichotomy. 

To prove Theorem 2.1 we must rule out the vanishing and dichotomy scenarios. 
Vanishing is ruled out as follows. Let denote a minimizing sequence. Then, 7i[ Uk ] = 
Z u + £fc, where e k \ as k — > oo. From the definition of H and the hypotheses Z„ < we have 

T v + e k = H[ u {k) } 

= ( -S 2 u« , u« ) - (a + l)- 1 J2 \^ ] \ 2a+2 



IrZ' 1 



i-l \,M\2a+2 
\ u l > 



and therefore 

\1„ | < v || 3« II,-. (7.6) 
Since vanishing implies || ||;°° — > as k — > oo, the lower bound ( |7.6| ) precludes vanishing. 
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Dichotomy is ruled out, as in the continuum case [22]|, using the strict subadditivity of the 
functional X v , i.e. if < a < v, then 

Z v < l a + l v „ a . (7.7) 
The idea is as follows. If dichotomy occurs (see (3) above) then as k — > oo and e — > we have 
J u = H[u {k) } + o(l) = H[a {k) ] + H[b^\ + o(l), (7.8) 
where we have used that aP^ and ¥- k > have disjoint supports. Furthermore, 

H[ a (fc) ] > I a+0 (e) and H[ 6 (fc) ] > l u - a+o{e) , (7.9) 
by definition of Ig Therefore, taking k — > oo and £->Owe get 

T v > l a + l v . a . (7.10) 

This contradicts ( [7.7|) . 
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